m 



o 



Finite amplitude inhomogeneous waves in 



(N : Mooney-Rivlin viscoelastic solids. 

^: 

^ ' Michel Destrade, Giuseppe Saccomandi 

<N : 2004 

o 

C/3 



Abstract 



B 

"i ■ New exact solutions are exhibited within the framework of finite 

p^ i viscoelasticity. More precisely, the solutions correspond to finite- 

O '_ amplitude, transverse, linearly-polarized, inhomogeneous motions su- 

perposed upon a finite homogeneous static deformation. The vis- 
coelastic body is composed of a Mooney-Rivlin viscoelastic solid, whose 
^ I constitutive equation consists in the sum of an elastic part (Mooney- 

OO ■ Rivlin hyperelastic model) and a viscous part (Newtonian viscous fluid 

!^ , model). The analysis shows that the results are similar to those ob- 

\^ I tained for the purely elastic case; inter alia, the normals to the planes 

f^ . orthogonal and conjugate with respect to the B-ellipsoid, where B is 

^^ ' the left Cauchy-Green strain tensor associated with the initial large 

static deformation. However, when the constitutive equation is spe- 
cialized either to the case of a neo-Hookean viscoelastic solid or to 
^ , the case of a Newtonian viscous fluid, a greater variety of solutions 

3 I arises, with no counterpart in the purely elastic case. These solutions 

include travelling inhomogeneous finite-amplitude damped waves and 
standing damped waves. 



1 Introduction 

Hayes and Saccoinandi [Ij initiated the study of finite amplitude waves in 
a simple class of viscoelastic materials of the differential type with a spe- 
cial constitutive equation for the Cauchy stress tensor. For these materials, 
denoted as Mooney-Rivlin viscoelastic materials, the Cauchy stress tensor 
is split into an elastic part, which coincides with the Cauchy stress tensor 
for a Mooney-Rivlin elastic material, and a dissipative part, which coincides 
with the Cauchy stress tensor for a Newtonian fluid (linear in the stretching 
tensor). Coleman and Noll [2] were the first to investigate the possibility of 
splitting the Cauchy stress into an elastic part and a dissipative part, and 
Fosdick and collaborators (Fosdick and Yu |.3J, Fosdick, Ketema, and Yu |1]) 
also considered Mooney-Rivlin viscoelastic materials. 

With respect to wave propagation, Hayes and Saccomandi |lj obtained all 
homogeneous plane waves that may propagate in Mooney-Rivlin viscoelas- 
tic materials maintained in a static state of pure homogeneous deformation. 
Then they showed [S] that principal transverse homogeneous waves may prop- 
agate when superimposed on a special class of pseudo-plane inhomogeneous 
steady motions. Finally, they considered antiplane shear waves superimposed 
on a static biaxiall stretch [6]. 

These results [H El E] provide a formidable corpus of exact solutions that 
may be useful as benchmarks for more complicated problems. Moreover, 
these solutions allow for a better understanding of the complex subject of 
nonlinear viscoelasticity ^ . The aim of this note is to extend these results 
to the case of inhomogeneous plane waves that is, plane waves for which 
the planes of constant phase are not parallel to the planes of constant am- 
plitude. Here we show that finite amplitude, linearly polarized, transverse, 
inhomogeneous plane waves superimposed on a generic homogeneous defor- 
mation are rare because two strict conditions must be met: first, the planes 
of constant phase must be orthogonal to the planes of constant amplitude; 
second, the equations of motion lead to a set of two differential equations 
to be satisfied simultaneously. However, when the constitutive equation is 
specialized from Mooney-Rivlin to neo-Hookean viscoelastic materials and 
further, to viscous Newtonian fluids, then these strict conditions disappear 
and many more exact solutions are obtained, including finite-amplitude in- 
homogeneous damped plane waves. A remarkable feature of the solutions is 
that although the static strains are finite, the viscoelasticity theory nonlin- 
ear, and the waves of a generic (separable) finite amplitude inhomogeneous 



travelling type, the governing equations nevertheless eventually reduce to a 
set of linear differential equations which are solved exactly. No small param- 
eters nor asymptotic expansions are required, in contrast with most studies of 
finite-amplitude waves in solids (see Norris [8] for a review). The only restric- 
tions lie with the specificity of the constitutive equation and with the form of 
the finite-amplitude inhomogeneous plane waves. This work follows the path 
laid out by Hadamard [S] for small-amplitude homogeneous plane waves in 
finitely deformed compressible elastic materials; by John pTO] and by Currie 
and Hayes [H] for finite-amplitude homogeneous waves in deformed elastic 
compressible "Hadamard" materials and incompressible Mooney-Rivlin ma- 
terials, respectively; and by Destrade [12] for finite-amplitude inhomogeneous 
waves in deformed elastic Mooney-Rivlin materials. 

The plan of the paper is the following. In the next Section we introduce 
the basic equations governing the constitutive model of a Mooney-Rivlin 
viscoelastic material, and the propagation of finite amplitude plane inhomo- 
geneous waves superposed on a large static pure homogeneous deformation. 
In Section 3, the equations of motion are solved and we find the conditions 
for the directions of polarization, of propagation, and of attenuation under 
which the waves may propagate. Sections 4 and 5 are devoted to the special 
cases of neo-Hookean viscoelastic materials and of Newtonian viscous fluids. 

2 Preliminaries 

2.1 Constitutive equations 

The elastic part of the model is characterized by the Mooney-Rivlin strain- 
energy density W which, measured per unit volume in the undeformed state, 
is given by 

2W = C{l~3) + D{ll-3), (2.1) 

where the constants C, D satisfy C > 0, D > or C > 0, D > ^^. The 
sum C + D is the infinitesimal shear modulus and I, II denote the first and 
second principal invariants of the left Cauchy-Green strain tensor B = FF : 

I = trB, 2H = (tr B)2-tr(B2). (2.2) 

The components of the gradient of deformation F are 



where Xi {i = 1,2,3) are the coordinates at time t of the point whose coordi- 
nates are Xj in the undeformed reference configuration. The Mooney-Rivhn 
viscoelastic sohd is incompressible, which means that only isochoric defor- 
mations are possible; in other words, the condition 

detF = l, (2.4) 

must be satisfied at all times. 

We assume that T^, the dissipative part of the stress, is given by 



^D 



z/(L + L^), L = FF-\ (2.5) 



where i/ > is a constant. The Mooney-Rivlin viscoelastic constitutive equa- 
tion for the Cauchy stress tensor T is given by 

T = -pl + CB-DB-^ + u{'L + l7), (2.6) 

where p is the indeterminate pressure introduced by the incompressibility 
constraint (12. 4p . to be determined from the equations of motion and eventual 
boundary conditions. 

We recall that in the absence of body forces, the equations of motion are 

div T = pi, (2.7) 

where p is the constant mass density of the material. 

2.2 Transverse inhomogeneous motions in a homoge- 
neously deformed material 

Consider a finite static isochoric homogeneous deformation defined by 

X = FX, (2.8) 

where the Fi^ are constant and det F = 1. For the deformation (12. 8p both 
B and B~^ are constant and L = 0. Then the following constant Cauchy 
stress tensor, 

T = -Pol + CB- DB-\ po = const., (2.9) 

clearly satisfies the equilibrium equations: div T = 0. 



On this state of static deformation, superpose a finite motion taking the 
particle at x to x, given by 

x = x + a/(b-x)^(n-x-t;t). (2.10) 

Here / and g are functions to be determined and n, b, and a are hnearly 
independent unit vectors, such that bxnT^O, b-a = n-a = 0. The motion 
f l2.10p represents a transverse inhomogeneous plane wave, propagating with 
speed V in the direction of n, linearly-polarized in the direction of a. When 
the wave is not damped, the amplitude varies in the direction of b; when the 
wave is damped, the amplitude varies in some direction in the (b, n)-plane, to 
be determined from the equations of motion. The case where / is a constant 
(homogeneous waves) has already been considered by Hayes and Saccomandi 
[1]; the case of an elastic Mooney-Rivlin material (z/ = in (12. 6p ) has been 
considered by Destrade [T2] . 

For the motion (I2.10p . the deformation gradient F is calculated using the 
chain rule as 

F = FF, F -.= 1 + f'ga®h + fg'a^n. (2.11) 

Taking the tensor product of F by 1 — fga ® b — fg'a n, we find the 
identity tensor, and so 

p-i=p-ip-i^ f-^ = 1- f'ga^h- fg'a^n. (2.12) 

It follows that the left Cauchy-Green tensor B and its inverse B are 
given by 

B = FBF^, B"^ = F-'^B-^F-\ (2.13) 

where F"""" is a shorthand notation for (F^^)^. Also, the time derivative of 
F is given by 

F = :^F, where F = -v{f'ga^h - fg'a® n). (2.14) 

Hence, the velocity gradient L is here 

L = fF~^ = fF^ = -v[f'g'a ® b + fg"a (g)n] = 't. (2.15) 

Now the corresponding Cauchy stress tensor T (say) takes the form, 

T = -pl + CB-DB'^ + i^(L + L^), (2.16) 

where the scalar p is to be determined from the equations of motion. 



2.3 Equations of motion 

In order to write the equations of motion, it proves practical to introduce the 
Piola-Kirchhoff stress tensor P associated with the motion fl2.10p and defined 
when the static state of pure homogeneous finite deformation is considered 
to be the reference configuration. Hence, 

Pi, := T,,F,-.i = -pF^^ + CF,,B,u - DBTj'Fj:/ + uP,,F^/ + uh^■ (2.17) 

and the equations of motion read 

pd^Xi/de = dP,k/dxk. (2.18) 

We introduce the real variables 

C:=b-x, 77:=n-x, so that / = /((), g = g(T]-vt). (2.19) 

Now we compute in turn the terms of (12.18^ . using the expansion ( 12.17^ 9 for 
Pik- First, the left hand-side term, 

pd%/dt^ = pv^fg"ai. (2.20) 

Next, the first term on the right hand-side, 

- d{pF^^)/dxk = -{dp/dxk)F^^ = -{p^ch +P,^nk)F^^ = -p^t^k -p^^rn, 

(2.21) 
where we used the Euler-Jacobi-Piola (EJP) identity for incompressible ma- 
terials dF^^ /dxk = in the first equality, and akUk = akbk = in the last. 
Now, the second term on the right hand-side, 

d{FipBpk)/dxk = {dFip/dxk)Bpk 

= [fgih ■ Bb) + 2/V(n ■ Bb) + fg"{n ■ Bn)]a,. (2.22) 

Then, the expansion of the third term on the right hand-side requires the 
EJP identity, 

d(B;j'F,/)/dxk = {dB.l/dx,)F,l = idB;l/dxk)6,, = dBT^/dx, 
= - [A(a ■ B-^b) + /V(n • B- V)]6, 
-[/V(a-B-ib)-f//(n.B-in)]n, 
-[f"g + 2fg\u-h) + fg"]Br^a, 

+ [2/7"/ + (3/'^ + fngg'ir, ■ h) + ff'{g'' + gg")]{s..B-'s.)h 

+ [2fg'g" + ffiSg' + gg")in • b) + (/'^ + fngg']ia-B''^)n,. 

(2.23) 



The expansion of the fourth term on the right hand-side also involves the 
EJP identity, 

= - v[f"g'hkhj + fg"{nkhj + hkUj) + fg"nknj\aiF^^^ 
= - v[f"g' + 2/'/(n • b) + fg"']a,. (2.24) 

Finally the fifth and last term on the right hand-side turns out to be zero: 

dF^^/dxk = -v[f"g'bkakbi + fg'nkakhi + fg'hkakrii + fg"'nkakni] = 0, 

(2.25) 
because akUk = akbk = 0. 

Now we proceed to the resolution of the equations of motion (I2.18P using 
( 1237]) and ([220!) -(ESS]). We treat in turn the cases of a Mooney-Rivlin 
viscoelastic material {C > 0,D > 0,1/ > 0), of a neo-Hookean viscoelastic 
material {C > 0,D = 0,1^ > 0), and of a Newtonian viscous fluid (C = 
0,D = 0,u>0). 



3 Mooney-Rivlin viscoelastic solid 

3.1 Orthogonality of n and b 

Destrade [12] proved that finite-amplitude inhomogeneous plane waves of 
evanescent sinusoidal type may propagate in a homogeneously deformed 
Mooney-Rivlin elastic material only when the planes of constant phase are 
orthogonal to the planes of constant amplitude. Here we prove that this 
result is not affected by the addition in (12. 6p of a viscous part to the Cauchy 
stress. 

Take the superposed motion a/(b ■ ^)g{n ■ x — vt) in (I2.10p to be of 
evanescent sinusoidal type, 

/ = ae-'^'"S g = 2/3 cos ujv-\ri-vt), (3.1) 

where the amplitudes a, (3, the frequency u, and the attenuation factor a 
are arbitrary real scalars. Then, leaving aside the pressure terms (I2.2ip for 
the time being, we find by inspection of fl2:20|) and fl2:22|) - fl2:25|) that all 
the other terms appearing in the equations of motion (I2.18P can be decom- 
posed along the set of linearly independent functions e"'^°'^ cos uv^^{ri — vt), 



e~^'^'' sin ujv^^{ri — vt), e^'^'^'^'^ cos2uiv^^{ri — vt), e~'^^'^'^ sm2uiv^^{ri — vt), and 
g-2wo-c fjj-^g only terms along this latter (time-independent) function come 
from the decomposition fl2.23p . where we have 

2f'f"g^ + ff'g' = -2a'pWai2a^ + v-^)e-^^^'^ + ... 

and ff'i3g'^ + gg") = -Aa'(5''uj''av-\-^''''^ + ... (3.2) 

Here the ellipses stand for terms proportional to time- dependent functions. 
It follows from the equations of motion that p must also be of this form, 

f = -a^/^^wVe-'""^ + . . . , (3.3) 

where pi is constant. Then the time-independent part of the equations of 
motion (I2.18P is simply 

= [-pi + D{2a^ + i;-2)(a-B'^a)]b + 2Dv-^{n ■ b)(a-B"^a)n. (3.4) 

Because n and b cannot be parallel for the motion (I2.10p to be inhomoge- 
neous, we conclude that pi = D{2a^ + f ~^)(a-B~^a) and that 

n-b = 0. (3.5) 

This important result, established in [12] for a Mooney-Rivlin elastic solid, 
still holds when the solid is viscoelastic with a constitutive equation of the 
form (12. 6p because as we just saw, the viscous terms (12.240 contribute to 
the equations of motion only with time-dependent functions. Turning back 
to a general inhomogeneous motion of the form (I2.10p . this result leads us 
to assume that here also, n and b are orthogonal, so that (n, b, a) is an 
orthonormal basis. Accordingly, some simplifications occur in the equations 
of motion. 

We introduce the notation 

pvl := C(n • En) + L)(a ■ B'^a) > 0, 
pvl := C(b • Bb) + D(a ■ B^^a) > 0, 

pb:=C{n-Bh). (3.6) 

Here f„ (vb) is the speed of an homogeneous finite-amplitude plane wave 
propagating in the direction of n (b) and polarized in the direction of a, as 



proved by Boulanger and Hayes [13]. Then the projection of the equation of 
motion f l2.18p along a is written in compact form as 

pvlf'g + 2phfg' - Mfg'" + f'g') + Pivl - v')fg" = 0. (3.7) 

Finally, the projections of the equations of motion f l2.18p along n and along 
b yield expressions for the derivatives p^ and p^; writing p^^ — P ri( gives 

m, 

if" 9 + f'9")i^ ■ B-^n) - ifg' + //')(a ■ B-^b) 

+ [(/'/" - ingg' - ig'g" - gg"')ff']i^ ■ B-^a) = o. (3.8) 

Hence the propagation of finite-amplitude, linearly polarized, transverse, 
plane waves in a deformed Mooney-Rivlin viscoelastic solid is governed by 
two equations: the "balance equation" (13. 7p and the "compatibility equa- 
tion" (13. 8p . to be solved simultaneously. 

3.2 Derivation of the solutions 

First divide (EZD by f{C,)g'{r] - vt) to get 

Pvlj^ + 2p6^ - uv{^ + ^) + p{vl - v'f-^ = 0. (3.9) 

Then differentiate first with respect to the argument of /, next with respect 
to the argument of g, to obtain 

{f"/f)'{g/g')' = 0, (3.10) 

so that either /"// or g/g' is a constant. The following cases cover all 
possibilities, 

{t)f"/f = -kl (h)/7/ = A;?, {^^^)g'/g = k2. (3.11) 

Without loss of generality, we pick Aji > in Case (i) and in Case (ii); 
in Case (in) we take A;2 > to avoid solutions which grow exponentially 
with time t (note that the forthcoming analysis can also be conducted to 
accomodate solutions which do grow exponentially with time). We treat 
each case in turn, keeping in mind that / must be trigonometric when g is 
exponential and vice-versa, for the motion to be inhomogeneous 

9 



Case (i): f{() = A cos ki( + B sin ki(. 
Substitution of this form of solution into the balance equation (I3.9p and 
differentiation with respect to ( yields ph{f' / f)' = and so, 

n-Bb = 0. (3.12) 

This condition means that the unit vectors n and b must be conjugate with 
respect to the central elliptical section of the B-ellipsoid, x ■ Bx = 1, by the 
plane orthogonal to a. Being orthogonal, they must be along the principal 
axes of this ellipse. 

Then the balance equation (13. 7p (with 6 = 0) gives a third order linear 
differential equation for g, 

^vg" - p{vl - v'^)g" - kfuvg' + kfpvlg = 0. (3.13) 

On the other hand, the compatibility equation (13. 8 p written as an identity 
for cosfciC, smki(, cos2A;i^, and sin2A;i^ yields 

g'g"-gg"' = 0, {g"-klg){si.B-'n) = 0, {g'" -kfg') {a-B-'h) = 0. (3.14) 

These three equations are satisfied simultaneously either when (a) g" = kfg 
or when (b) g"/g =const., a ■ B~^n = a ■ B~^b = 0. 

In Case (ia), g is an exponential function. Discarding solutions which 
blow up with time, we find 

g{ri-vt)=e'''^''-'"\ (3.15) 

Substitution into the third-order differential equation (I3.13P fixes v as 

pv' = p{vl - vl) = C[(n ■ Bn) - (b ■ Bb)]. (3.16) 

The quantity v is real when n and b are in the respective directions of the 
minor and major axes of the elliptical section of the B-ellipsoid by the plane 
orthogonal to a. 

In Case (ib), g is of the form 

g(^ri-vt)=e^'^''-^'\ (3.17) 

where ^2 > to avoid blowing-up solutions. The conditions on n, b, a, 
together with (I3.12p . imply that these unit vectors are along principal direc- 
tions of the B-ellipsoid. Substitution into the third-order differential equa- 
tion (I3.13P relates ki to k2 through 



pvl - k2VV 



kl = "^^"..2i;..r^ ^2- (3.18) 

b 
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Recall that ki is assumed real so that, in order to construct a solution here, 
we may pick any values for k2 > and for w > 0, as long as p{v^ — v^) — k2i'v 
and pv^ — k2i^v are of the same sign, where f„, Vb are given by (13 .Gp with n, 
b, a along principal directions of the B-ellipsoid. Then ki is given by fl3.18p . 

Case (a): f{() = A cosh kiC. + B sinh kiC^. 
Substitution of this form of solution into the balance equation fl3.9p and 
differentiation with respect to C yields the condition fl3.12p as in Case (i). 

Then the balance equation f l3.7p (with 6 = 0) gives a third order linear 
differential equation for g, 

^vg" - p{vl - v^)g + k\vvg - k\pvlg = 0. (3.19) 

Also, the compatibility equation (13. 8p written as an identity for coshA;i(^, 
sinhfci^, cosh2fci(^, and sinh2A;i(^ yields 

g'g"-gg"' = 0, ig"+klg){^.B''n) = 0, {g"'+klg'){si.B-'h) = 0. (3.20) 

These three equations are satisfied simultaneously either when (a) g" = —kfg 
or when (b) g"/g =const., a ■ B^^n = a ■ B^^b = 0. 
In Case (iia), g is of the form 

g{ri — vt) = di cos kiij] — vt) + d2 sin kiij] — vt), (3-21) 

and substitution into the third-order differential equation (I3.19P fixes v again 
as ( KT^ .v^ = vl-vl 

In Case (iib), g is of the form g = diCOsk2{r] — vt) + d2smk2{i] — vt), 
where /c2 is arbitrary. The conditions on n, b, a, together with (I3.12p imply 
that these unit vectors are along principal directions. Then substitution of g 
into the third-order differential equation (I3.19P gives 

uv{kl - kl)g' - p[kX - kl{vl - v^)]g = 0. (3.22) 

This equation is satisfied for g trigonometric and v ^ Q only when the respec- 
tive coefficients of g and g' are zero, conditions which lead back to a solution 
of the form (I3.2ip and v uniquely determined by (I3.16p . The existence of 
these 'special principal motions' is therefore more limited than when the 
solid is purely elastic ^^ because there they may propagate at an arbitrary 
speed V within the interval [0,f„]. 

Case (m): g{ri - vt) = e^^^''""*). 
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Substituting this form of solution into the balance equation fl3.7p . we 
obtain a second-order linear differential equation for /, 

{pvl - huv)f" + 2hphf + kl[p{vl - v"") - huv]f = 0. (3.23) 

Also, the compatibility equation fl3.8p reduces to 

if" + klf){^ . B-^n) - Uf" + klf){^ . B-^b) 

+ k2{ff" - ff"'){a ■ B-^a)e''^^'^-'"'> = 0. (3.24) 

Differentiating this equation with respect to (r/ — vt), we find that /"// = 
const. This constant must be negative so that the combination of / and g 
represents an inhomogeneous motion [12]. But if / is trigonometric, then 
the second order differential equation (I3.23P can be satisfied only when pb = 
C(n-Bb) = 0, which leads us back to Case (i) and the solutions (I3.15p - (l3.16p 
(non-principal motions) and (I3.17p - (l3.18p (principal motions). 

3.3 Summary of results and comparison with the purely 
elastic case 

For a deformed elastic Mooney-Rivlin material, Destrade [I2] proved that 
only two types of finite-amplitude inhomogeneous motions are possible: ei- 
ther / is trigonometric and g is exponential, 

/(() = ai cos kiC + 02 sin kiC, 9{v - vt) = rfie^^^""^*) , (3.25) 

or / is hyperbolic and g is trigonometric, 

/(C) = fli cosh ki(+a2 sinh ki(, g{rj—vt) = di cos k2{rj—vt)+d2 sin k2{rj—vt), 

(3.26) 
where ai,a2,di, ^2 are constants and /c2 is arbitrary. 

When the vectors (a, b, n) are not aligned with the principal axes of 
deformation, then the quantities ki and v are determined by 

ki = k2, pv^ = p{vl - vl) = C[{n ■ Bn) - (h ■ Bh)], (3.27) 

and the unit vectors n and b must be conjugate with respect to the B- 
ellipsoid, x ■ Bx = 1, that is the condition (I3.12p must be satisfied. 
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When the vectors (a, b, n) are ahgned with the principal axes of defor- 
mation ('special principal motions'), then the quantity v may take any value 
within the interval [0, f „] and ki is given by 

kl = "^^^^kl (3.28) 

Here we saw that the corresponding results for a deformed viscoelastic 
Mooney-Rivlin material are essentially the same as above, with the following 
difference. When the three vectors a, b, n are not aligned with the principal 
axes of deformation, then the solution is indeed of the form (I3.25p - (l3.27p with 
the condition (13.121) . However, when all three vectors a, b, n are aligned with 
the principal axes of deformation, then the solution is of the form (13.251) or 
(I3.26P but: when / is trigonometric and g is exponential as in (I3.25p . then 
kl is given by (I3.18p . and v, ^2 are arbitrary as long as kf is positive; when 
/ is hyperbolic and g is trigonometric as in (13.261) . then v are determined by 
(I3.16P that is, v is fixed and no longer arbitrary as in the elastic case. 

Now we see that when D = 0, the visco-elastic solid allows more solutions 
than its elastic counterpart. 

4 Neo-Hookean viscoelastic solid 

When D = {C y^ 0) the strain-energy density (12. ip reduces to the well- 
known neo-Hookean form. A great variety of exact solutions is uncovered, in 
particular solutions in the form of finite-amplitude damped inhomogeneous 
waves. Note that this problem was recently studied in the elastic compress- 
ible case by Rodrigues Ferreira and Boulanger [H]. 

We start with the equations of motion (I2.18P written at D = 0. The 
unit vectors b and n are not necessarily orthogonal but the triad (a, b, n) is 
composed of linearly independent vectors, so that the coefficients along each 
vector in the equations of motion are all zero. The coefficients along b and 
n are respectively, 

= -p^^, = -p,^, (4.1) 

from which we deduce that p =const. satisfies this part of the equations of 
motion. The other part, along a, is independent of p and reads 

p{v' - vDfg" - pvlf'g - 2pbf'g' + uv[f"g' + 2fg"{n ■ b) + //'] = 0, (4.2) 

13 



where the quantities 

pvl = Cn ■ Bn > 0, pv^ = Cb ■ Bb > 0, pb = Cn- Bb, (4.3) 

were used. 

Now divide f l4.2p by f{C)g'{i] — vt) and differentiate, first with respect to 
C and then with respect to {rj — vt), to obtain the necessary condition 

We separate the functions of different variables in this equation to conclude 
that either 

The following cases cover all possibilities, 

(^) /' = hf, (2z) g' = hg, (4.6) 

and (in) pv^f" — kif — k^f = 0, 2z/t>(n ■ h)g" — k2g' — kig = 0, where ki, ^2, k^ 
are constants. In fact, using expressions for /", g", and g"' derived from {Hi), 
and substituting into (14.21) . it is straightforward to check that Case {iii) leads 
to Cases {i) or (ii). 
In Case (i) we have 

/(C) = exp(fciC), 

g'" + [A(^2 _ ^2) ^ 2fc^(n • b)]/ + fci(fci - 2-^b)g' - p^klg = 0. (4.7) 
vv vv uv 

In Case {ii) we have 

gijq - vt) = exp[k2{r] - vt)], 

(pvl - k,uv)f" + 2k2[pb - uvk^in ■ h)]f - kl[piv^ - vl) + vvk2]f = 0. 

(4.8) 

Now consider (14. 7p in greater detail. The characteristic equation associ- 
ated with the third-order differential equation for ^f is a cubic. Depending on 
the values of the coefficients (that is, on the choices for v and ki), the cubic 
has either only real roots or one real root and two complex conjugate roots. 
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Real roots lead to an homogeneous motion and must be discarded. Pure 
imaginary roots lead to inhomogeneous plane waves with attenuation in the 
direction of b and propagation in the direction of n. Non-pure-imaginary 
complex roots lead to inhomogeneous damped plane waves with attenuation 
in the plane of b and n, exponential damping with time, and propagation in 
the direction of n. Because v or ki may be prescribed a priori, there is an 
infinity of such solutions, of which we present one explicitly. 

We choose v such that the coefficient of g" in fl4.7p 9 is zero, that is we fix 

V as 

pv = pvo := kiu{n ■ b) + ^[kiu{n ■ b)]2 + (pw„)2. (4.9) 

Then (14. 7p reduces to 

9"' + (3o9'-lo9 = 0, ^o = h{k,-2-^), ^o = kl^. (4.10) 

l/Vo Wo 

Now we seek a solution in the form 

g[r] - Vot) = e^°(''-'^°*) cosWo(r/ - v^t), (4.11) 

where Aq and Uo are real constants determined as follows. Substituting (14.111) 
into (I4.10p . we find in turn that Uo and Xo satisfy 

Lol = 3\l + /3o, 8\l + 2(3oXo + 7o = 0. (4.12) 

The real root of the cubic in Xo is 

Xo = 3^ [12^12/33 + 8l7o-1087o]^ -/3o[12v/12/33 + 8l7o-1087o]"i (4.13) 

That (jJo is real is ensured for instance when /3o > (see Eq. (I4.12p i). A few 
lines of calculation show that a sufficient condition for Po to be positive is 

,?> mi = ApCin.Bhf 

' z/2[t;2+46(n-b)] z/2[(n-Bn)+4(n-Bb)(n-b)]' ^' ' 

Thus for any arbitrary value of fci satisfying this inequality, the following 
finite-amplitude damped inhomogeneous plane wave may propagate in a de- 
formed neo-Hookean viscoelastic body, 

X = X + aae^k^h+Xony^-XoVot cosa;„(n • x - Vot), (4.15) 
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where x is given by f l2.8p . a is arbitrary, the orientation of the triad of unit 
vectors a, b, n (b and n orthogonal to a) is arbitrary, and Aq, Wq, and Vq are 
given by (14.131) . (14. 121) i . and (14. 9p . respectively. The wave propagates in the 
direction of n, is exponentially damped with time, and is attenuated in the 
direction of fcib + Aon. 

To conclude this Section, we point out that inspection of (14. 4 p makes it 
clear that a neo-Hookean viscoelastic constitutive equation yields more in- 
homogeneous waves than a neo-Hookean elastic constitutive equation, where 



5 Newtonian viscous fluids 

With C = D = 0, equation (12. 6p is the constitutive relation defining an 
incompressible Newtonian viscous fluid of viscosity v^ and (12. 7p . or their 
specialization (14. 2 p at C = 0, are the (Lagrangian) Navier-Stokes equations 
in the absence of body forces for the finite amplitude inhomogeneous plane 
wave (12301) . that is fg' + [2(n ■ b)/' + {pvlv)f\g" + fg'" = 0. Integrating 
with respect to the argument of g and taking the constant of integration to 
be zero for simplicity, we obtain 

f"g + [2in-h)f' + ^f]g' + fg" = 0. (5.1) 

With methods similar to those presented in the two previous sections, it 
is a straightforward matter to show that either f or g are pure exponential 
functions that is, either 

/(C) = exp(A;iC), and g" + [2k,{n-h) + ^]g' + k^g = 0, (5.2) 

or 
g{7^-vt) = exp[k2iv-vt)], and /"+2A;2(n ■ b)f+A;2( — +^2)/ = 0. (5.3) 

where ki, ^2 are arbitrary constants. The latter case represents a standing 
damped wave, whatever the form of / may be. The former is a travelling in- 
homogeneous damped wave as long as the characteristic equation associated 
with the second-order differential equation for g has complex roots. This 
condition reads 

< pv < 2kiv[l - (n ■ b)] =: pvo. (5.4) 
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Hence the constants ki and v are arbitrary as long as ki is positive and v 
is less than vq defined above. The following displacement field u is solution 
to the Navier-Stokes equations for any orientation of the pair of unit vectors 
b, n in the plane a ■ x = (which may be taken without loss of generality as 
z = say, because the fiuid is isotropic), 

Ml = 0, U2 = 0, 

where A, B are constants and k is defined by 



l-[(n.b) + -^p. (5.6) 

This finite-amplitude wave propagates in the direction of n with speed v 
and wave number /cik, is damped exponentially with respect to time, and 
is attenuated in the direction of b + a/1 — ft^^n. The angle 9 (say) between 
the normal to the planes of constant phase and the normal to the planes of 
constant amplitude is given by 

cot = n ■ b + Vl - fi;2 = 2(n ■ b) + -^. (5.7) 

kiu 

Note that Boulanger, Hayes, and Rajagopal [15] also investigated the pos- 
sibility of linearly-polarized, damped, inhomogeneous plane wave solutions 
to the Navier-Stokes equations. However, they prescribed a priori the form 
of the solution. The result presented above allows for a greater variety of 
solutions, in the same manner that starting from a finite-amplitude wave 
in the form (I2.10p where / and g are unknown functions, rather than pre- 
scribed as exponential and sinusoidal functions, allowed Rodrigues Ferreira 
and Boulanger [1^ to find more solutions than Destrade [16] in the elastic 
neo-Hookean case. 

We also note that the solutions found here for the Navier-Stokes equa- 
tions are a generalization of the celebrated Kelvin modes [T^. These are 
disturbances in the form of homogeneous plane waves of particular interest 
in the stability study of basic fiows characterized by spatially uniform shear- 
ing rates, i.e. self-equilibrated homogeneous motions [18]. This is because 
Craik and Criminale [19] have shown that the superimposition of a single 
Kelvin mode to the above mentioned basic fiows is an exact solution of the 
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full Navier-Stokes equations. It would be interesting to investigate if the 
superimposition of our solutions to homogeneous self-equilibrated motions 
will give some new exact solutions for the Navier-Stokes equations, but this 
investigation is outside the framework of the present paper. 
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